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DIFFERENTIAALVERGELIJKINGEN 


Grondbegrippen 


Tenzij anders vermeld, zijn in het vervolg alle variabelen 

en functies reel. 

Zij D ei & open intervallen van de reêle rechte, F een gege- 
ven functie (x,y) * Flxsy) » Xx € D, y es. 

Zij voorts V de verzameling continu dirtenentieenbace functies 
x * y(x), x E Díd.w.z, V = Bt. Wij fennie het volgende probleem: 


Bepaal een element y € V zodanig dat 


xe) UN - p(x‚yo). 


De uitdrukking 


dy - 
EE F(x,y) 


noemen wij een differentiaalvergelijking (van de eerste orde). 

Een functie y(x) die voldoet aan de eisen van ons probleem noemen 
wij een oplossing van de differentiaalvergelijking. 
Probleemstelling zoals boven gegeven kan op verschillende manieren: 
worden gegeneraliseerd. 

Enerzijds, kunnen wij n-de orde differentiaalvergelijkingen be- 
studeren. Er wordt dan gevraagd functies y(x) € P(D) te bepalen 


die voldoen aan 


n n= 
Sd =F kre, EEE: DN ad _9 
dx dx dx 


waarbij F wederom een gegeven functie is. 


Anderzijds kunnen wij stelselsvan eerste orde vergelijkingen 


onderzoeken. Men vraagt dan naar stelsels van functies 


yy syg lx), ee y(x) die voldoen aan 





de E_(xsy4- . Yn) 


In dit college zullen wij ons hoofdzakelijk beperken tot n = 1 


f - ) 
| Ee = Foy gore syn) 

| dy, 

| ran FE (xsyYqores Yo) 

| dy, 


en n= 2. 
Er bestaat een nauw verband tussen hogere orde diff.vergelijkingen 
en stelsels van eerste orde vergelijkingen. Dit kunnen wij 


Inzien als volgt. 


Beschouw 
7, 
dy = Poe 
2 dx 
dx 
Definieer 





- „BE 


Er ontstaat dan een stelsel 





Ë dy, 


| | dx” /2 


dy, 
oe B 





Wij noemen het stelsel eerste orde diff.vergelijkingen lineair, indien 


en ee ee oh HS De weken a 
NS 


u M5 


Foy qse 9) = aj (x-y3 + b;(x) bn SS Meens 


J=1 
Het stelsel is niet-lineair als aan de eisen van lineariteit 


niet is voldaan. Een n-de orde differentiaalvergelijking is 








is lineair als het corresponderende stelsel eerste orde differenti- 


aalvergelijkingen lineair is. 


Het beginwaarde-probleem 


In toepassingen zijn de functies y, of Vg" °Yj vaak te inter- 


n 
preteren als grootheden die de toestand van een systeem beschrijven, 
terwijl in de differentiaalvergelijkingen de wetten zijn ge- 
formuleerd volgens welke het systeem functioneert. Veelal is 


de toestand voor zeg Xx = Xx, bekend en men wenst te bepalen hoe 


0 


de toestand voor x # Xx, zal zijn. Er ontstaat aldus het zoge- 


naamde beginwaarde-probleem : 


n_= 1 Bepaal y(x) zodanig dat 


0) 
Se z F(xoy) : ye) = y 


n = 2 A.Bepaal yy sy (ax) zodanig dat 


de Ng 


es ed 
yo’ ie yy ; Yo xp) = Y 5 


ax 5 F‚lsy4sy9) 


n = 2 B.Bepaal y(x) zodanig dat 


an 


) | 
dy z r(x,ys) 5 ylang) = yo 3 (42) = yy 
ee dd P 0 2 dx En 2 


Het is à priori niet zeker dat een beginwaarde-probleem een 


oplossing bezit. Beschouw als voorbeeld 


ae 


1 
X 


zi USS Es vl == Le 
Het is eenvoudig in te zien dat alle functies die aan de 


differentiaalvergelijkingen voldoen gegeven zijn door 


y(x) = Inx + c 





maar in deze familie van functies is er geen die voldoet aan 

y(0) = 1. 

Anderzijds, zo er oplossingen bestaan, moet men de vraag stellen 
hoeveel er bestaan. Uit het oogpunt van toepassingen op deter- 
ministische systemen zou men wensen dat er slechts Één oplossing 
bestaat. Maar niet alle differentiaalvergelijkingen vervullen 
deze eis. 


Beschouw bijvoorbeeld 


= Wy; 0Sx<SL<1 3 y(O) = 0 


ed 


Er voldoen twee oplossingen, te weten 


y(x) = 0 en y(x) = H. 


Wij formuleren nu, voor n = 1, een algemene existentie en een- 


duidigheid stelling. 


Definitie. De functie F(x‚y), x ED ; y EG, heet Lipschitz 
continu t.a.v. y in D XG als er een constante À 


bestaat zodanig dat 


(1) 


(Y xE DD) (Y y (2) 


E G)(Y y' TLE G) 


rx) — rixyt 2) satytt? - yer. 


Stelling. Zij Flx‚y) continu t.a.v. x en Lipschitz continu 


Erdee Van D X 5. 


Voor elke yv E G bezit het probleem 


Û 


d - | - 
dx F(x,y) ì vx) y 








Eén en slechts één oplossing y(x), gedefinieerd 


in een interval |x - xl S bels > 


Opmerkingen. Bewijs van deze stelling zullen wij in een later 
paragraaf schetsen. De stelling is gemakkelijk uit te breiden 
tot stelsel met n 2 2, er moeten’ dan aan elke F.si tem 


analoge eisen worden gesteld. Wij verwijzen hiertoe naar het 


college "grondslagen van differentiaalvergelijkingen". 


Enkele elementair integreerbare vergelijkingen. 


Bij de studie van differentiaalvergelijkingen is een aantal 
elementair oplosbare classen niet alleen een zinvol uitgangs- 
punt, maar ook een bijzonder nuttig hulpmiddel bij verdere 


onderzoekingen. 


II. De autonome vergelijking van de eerste orde. 


Een vergelijking van eerste orde heet autonoom als zij de 


gedaante heeft 


end = g(y) 


ZAj nu Yo zodanig dat EY) F 0. Wij definiëren 





ki 

Jo 
d _ 1 dy (x) 
rde bdhadd ad 5772) Cam 


Als y(x) oplossing is van de differentiaalvergelijking, dan geldt 


d B 
oe G(y(x)) = 1 


waaruit volgt Gly(x)) = Xx + Cc 
met c een willekeurige constante. 


Wij hebben aldus een relatie gevonden ‘tussen y en x. Hieruit 

















moet alsnog y(x) expliciet worden bepaald. 


dy . 
Voorbeeld : Ees y 
G(y) = In 4- 5 In AIEE X + C 
n Yo Jo 
Als y(x) SV dan is c = -Xg° Invoeren van inverse functie 


® 


levert uiteindelijk 
(x-xg) 
y(x) = Yp © | 


II. De vergelijking met gescheiden variabelen. 


Wij beschouwen nu 


dy - 
Ee f(x) g(y) 
X 
Za E == í £Cxijdxs 
X0 


Er volgt onmiddelijk 


TE = g(y) 


De vergelijking kan nu worden opgelost als onder I. De alge- 
mene relatie tussen y en x 1s 
dn Xx 


AKE el Sasi * GQ 
Ô ê 


5 


Door deze relatie is y impliciet als functie van x gegeven in 
een omgeving van Xp 


(Ga na): 


III. De lineaire inhomogene vergelijking 


De lineaire vergelijking 


nd = f(x).y + ò(x) 


noemt men inhomogeen als ò # 0. De vergelijking is niet van 


het type van gescheiden variabelen. 





Wij passen toe de volgende transformatie (genaamd “variatie 
van constante!) 


y= nlx).0G0) 
met n(x) een oplossing van het homogene probleem 


an 


Ee f(x).n 


Wij kiezen 


E ded 
n(x) = exp {f f(x )dx} 
X 


0 
Invullen levert 
de _ 
n(x) ers dx) 
X 1 eZ 
En bijgevolg 0x) = f Oe dxt + Cc. 
Xo nx) 


IV Lineaire homogene stelsels met constante coefficieëenten. 


dy ; n n 
8 er nnn ‚> “.s e 1 = nde dà 
Beschouw Er oe a;; y; 1 ds sn 
Jed 
waarbij a... constanten zijn. Als n = 1 dan zijn de oplossingen 


1 
exponentiële functies. Het lijkt de moeite waard te onderzoeken 
of exponentiële functies ook oplossingen kunnen zijn voor n > 1. 
Stel daarom y; (x) = Bn LS satans MH 


Er zou dan moeten gelden 


n 
X a;ses T Ae; Ek Ok Aang Ù 


Dit stelsel lineaire algebraische vergelijkingen voor c,,..: +. 


bezit niet triviale oplossingen als het coêficiënten determinant 


gelijk nul is. 


Voor n = 2 vinden wij aldus 
aa Te 
= 0 
nl Sh 4 


DT en ee en etste amen ve ven meereden va er dee 





EERENS DE UAR OO OE OR MR 





en enn ne ie er de nn ne arr eer ee kn a 





Dit is de eigenwaardevergelijking behorende bij de coëfficiënten 
matrix eed De wortels van de resulterende tweedegraads 


vergelijking voor \ bepalen de gezochte exponentiële functies. 


Locaal onderzoek : iteratie en reeksontwikkeling 


Beschouw nu; voor n = 1, het algemene probleem 


zie 


Ely) 3 yCxg) = Yo 


Het probleem kan equivelent worden geformuleerd als de inte- 
graalvergelijking 
8 | 
y(x) = Yo * S_ F[E‚y(E)ldE 
pd 
0 
Volgens het iteratie proces (van Picard) definiëren wij een 
rij functies 
5 | 
= | ee 
yo) Yo Ek F[E‚y__4 (E18 ne 4 
X 
0 
Het proces is constructief, d.w.z. elke y(x) kan worden be- 
rekend. De hoop en de verwachting is dat y(x) convergeert 
(als n > @) en wel naar de functie y(x). 
De rij y kan (op een nagenoeg triviale manier) worden 
geïdentificeerd met 
n | 
z _ ° > 
Ie Ig t E yo) Yp} e Med 
p=1 
Wij onderzoeken nu de convergentie van de reeks, die ontstaat 
als n > @, 


Beschouw (om de gedachten te bepalen) x > Xr 


X 
ly, - Yp-1C9| SS |FLEsy_4 (8) + P[E,y,_(E)llde,p > 2. 


X0 








Wij veronderstellen 
F(x,;y) continu in D XG, en in het bijzonder Lipschitz 
contanu tsasVa Ye 


Xx 


En ef R 
0 D 45 Yo G 


Wij merken op dat, voor alle n, y(x) continue functies zijn 


met ACR, ne Bijgevolg bestaat er een interval x, <x $ Xx» 
X4 > Xp» zodanig dat y(x) E G. Toepassing van Lipschitz 
continuïteit levert in Xo S X SX, 3 


5 | 
y,6) u Yp-1 91 SA f yj-168) — y (E)IdE 5 p= 2. 
x 


p-2 
Q 


Merk nog op dat 


x 
yv, Go E Yo "4 F(Esy garde. S M(x-x 
8 | 


waarbij M een constante. 
Met volledige inductie vindt men: 


N 5 k EN w P 
<M [Ax xo’] 
NE Ne 


Cx) | Se 


yo oe Jp-1 


Er volgt dat de termen van de reeks 
‚ es 
oe ly) - Yp-1C9} 
in absolute waarde term voor term worden gemajoreerd door 
termen van een convergente reeks. Bijgevolg convergeert ye) 


en wel gelijkmatig in elk interval Xo S Xx S Xx, waarvoor geldt 


1 


dat gn (xx) S CG. Er kan nu nog worden aangetoond dat inderdaad 


lim y, aan de differentiaalvergelijking voldoet. 
noo 


Aldus is constructief existentie van oplossingen aangetoond. 
Zij tu en, = lim ve 9). om de eenduidigheid van ons probleem 
noo 


te onderzoeken shellen wij da 5 Ge) een ändere functie is die 


voldoet. 








Met afschatting analoog aan het voorgaande vinden wij 
5 | | 
y,G) gl SAS [Yr e BCEIE A med 
ze | 
Q 


Ly (0) - To) S M(x-x,) 


EL VOLGE 


M [Alx-x ltd 
Ly, Go) ad T(x)| S X TCn+ti) E 


Dus: lim ly Go - EGO] = y(x) = Tx) = 0, 
n 


=> OO 


d.w.z. E(x) is identiek met vx). 

Het bovenstaande iteratie proces, en bewijs van convergentie, 
kan worden uitgebreid tot stelsels van willekeurig vele eerste 
orde differentiaalvergelijkingen. Aldus is zowel een theoretische 
basis voor de beginwaarde problemen gelegd, als een practische 
mogelijkheid geschapen om de oplossingen te construeren. Echter, 
in toepassingen zijn de practische mogelijkheden van het construc- 
tief proces nogal beperkt. | 

In het algemeen is het berekenen van y, Voor grote waarden van 

n zeer moeizaam en arbeidsintensief, terwijl de uitkomsten 
onoverzichtelijk zijn, zodat men gedragingen van y(x) nauwe- 
lijks nog kan onderzoeken. | 
Practisch gesproken is het constructief proces in het algemeen 
slechts bruikbaar voor locaal onderzoek, d.w.z. | x-x, | klein. 

Dan convergeert y, snel, en met n niet al te groot kan al een 
redelijke benadering van y worden gegeven. 

Het analyseren van de oplossingen van differentiaalvergelijkingen 
in grotere intervallen XS x Sd (in gevallen die niet ele- 
mentair oplosbaar zijn) blijkt dan ook een voortdurende uit- 


daging voor de toegepast wiskundige. 


Autonome stelsels. Baankrommen in het fasevlak. 


Voor een classe van de zogenaamde autonome differentiaalverge- 
lijkingen is het mogelijk met speciale methoden het gedrag van 
de oplossingen te analyseren. 


Wij beschouwen het stelsel 


Lr EP PREM Yo) 
dx Ne nd NAE 2 Vd ad 

f‚ en f), voldoen aan de eisen van de existentie- en eenduidig- 
heid stelling. 

Zij y(x): y(x) een paar functies die voldoen. 

Men verifieert gemakkelijk het volgende 
As yi): y(x) voldoet, dan voldoet ook y4(E) sy, (ED), met 

5E =x= Xp Xo willekeurig (maar vast). | 

Wij gaan nu het gedrag van de oplossingen van het stelsel 
onderzoeken in het Y47Y vlak, dat wij het fasevlak noemen. 
Elk paar niet constante functies y‚(x),y,(x) die aan de diffe- 
rentiaalvergelijking voldoet definieert in het fasevlak een 
verzameling punten die wij een baankromme noemen. 

De existentie- en eenduidigheid-stelling voor oplossingen van 


beginwaarde problemen laat zich vertalen als volgt 


Stelling : Door elk punt Y4°Yj E G gaat één en slechts één 
baankromme. 
Er kunnen in G punten voorkomen met speciale en uitzonderlijke 


eigenschappen. 
Definitie: (n‚‚n,) heet een singulier punt (of estetsk punt, 

es, stationnair punt) als £ ‚m4 on) = Ë (ns on) = 0, 
Eigenschap: Zij (nson,) een singulier -punt (in G). 

Als yo) = Ms le) = n dan geldt voor alle X 


ON Sena Hele) So 














RE ee NEEDE NEP Te OTE ET RNB MPEG NRR, 


Opmerking : Een singulier punt is dus een baankromme die “tot 

Êén punt is gedegenereerd. 
Buiten de singuliere punten kan een differentiaalvergelijking 
voor de baankrommen worden opgesteld. 
Het is eenvoudig in te zien dat geldt 

Ay, _ falFgsYa) 

Onderzoek van baankrommen in het fasevlak heeft twee belang- 
rijke aspecten. Enerzijds kan soms de differentiaalvergelij- 
kingvan de baankrommen tot een gemakkelijk integreerbaar type 


behoren. Anderzijds, in het meest algemene geval, kan daar zo- 


genaamd kwalitatief onderzoek, vergaand inzicht worden verkregen 


in het gedrag der oplossingen. Wij bespreken hier het tweede 


aspect. 

De differentiaalvergelijking der baankrommen definieert in het 
fasevlak een richtingsveld, d.w.z. : in elk punt Y4°Y7 is de 
helling van de baankromme die door dit punt gaat gegeven. 

De richting waarin een punt yy) lands een baankromme 
“beweegt! voor bijv. toenemende x kan nog uit de oorspronkelijke 
differentiaalvergelijking van het stelsel worden afgeleid. 
Men kan nu proberen bij benadering de baankrommen te schetsen 
door in het richtingsveld "zo goed mogelijk passende!" krommen 
te CONSTruUeren. | 

Als een nagenoeg triviale illustratie van bovan beschreven 


procedure beschouw 


Vi U 
dx 1 * dx 72 
Er is één singulier punt: Y, * Y, * 0. De vergelijking van de 


baankrommen is 





dy, Yy 


dy, Y4 


Richtingsveld is in onderstaande schets aangegeven. 


2 





Het is evident dat "zo goed mogelijk passeerde" krommen rechte 

„lijnen zijn; dit zijn in dit geval ook de exacte oplossingen 

van de vergelijking voor baankrommen. 

Beschouw thans, als een meer serieuse toepassing, het stelsel 
dy 4 dy 
A he eN 


Er zijn twee singuliere punten, te weten 
Me Oh Ae 
De vergelijking van de baankrommen is 


dy y,{1-y,) 

dy, y4ldey) 
Men ga na dat het kwalitatief onderzoek een "“faseportret! van 
de oplossingen oplevert, zoals in nu volgende schets gegeven 


(pijltjes geven aan "beweging! voor toenemende x). 








hindi beh aante neee ien area re erken as die bad pr aant and aus haken rn meoctsr Seher 
rd % 8 





Ya 





Wij merken op dat de vergelijking van de baankrommen van het 
type van gescheiden variabelen is. Men kan dus, volgens het 
eerder besprokene, integreren, maar uit het resultaat is in 


dit geval bijzonder moeilijk Y) expliciet in 4 uit te drukken. 


De lineaire vergelijking van tweede orde. 





We beschouwen nu de lineaire vergelijking 


2 
dy dy = 
de + b(x) en e(x)y 0 


Op een wijze die eerder was aangegeven kunnen wij met deze 


vergelijking associeren het stelsel 


dy dy 
1 REE 3 
rm D(x)y, c(x)y, 





ee 


Zij nu b(x) en e(x) continu en begrensd voor x € D. Dan zijn 

de rechterleden van het stelsel Lipschitz-continu t.a.v. Y4 

en yY, Op de gehele reêle rechte, en aan de eisen van existentie- 
en eenduidigheid-stelling voor beginwaarde problemen is voldaan. - 
Wij definiëren twee speciale oplossingen ò(x) en W(x) voor de 


differentiaalvergelijking als volgt 


2 
d dy n e de 5 
f + De + cb = 0 ò Cx) s Js (st) _ = 0 
8 | . 
dp dp n : du 5 
dx ‘ XX | 


Een willekeurig beginwaarde probleem wordt nu opgelost door 


y(x) = A Wlx) + B plx) 


waarbij vx) A 3 (de ) _ B. 
X=X 


Anderzijds, met A en B willekeurig, stelt de uitdrukking hier- 
boven de verzameling van alle functies voor die aan de differenti- 


aalvergelijking voldoen. Een dergelijke uitdrukking noemt men de 


» 


algemene oplossing. 


Opmerking: Dat er inderdaad geen andere functies bestaan die aan 
de differentiaalvergelijking voldoen wordt aangetoond door de 


nu volgende overweging: stel een functie u(x) voldoet aan de 


differentiaalvergelijking. Men kan dan uitrekenen u(x) z A 3 
(2) = B. Echter de functie. 
X=X 


adlx) + BWCXxX) 
voldoet eveneens aan de ieenttaawesderidkane en aan de- 
zelfde beginwaarden. Wegens eenduidigheid volgt 

u(x) = adlx) + Bx) 


d.w.z. u(x) is bevat in de algemene oplossing. 





_ 16 =— 


Het randwaardeprobleem. 


Wij formuleren nu een nieuw probleem voor de tweede-orde 
differentiaalvergelijking 


Bepaal een oplossing y(x) in Xp S x S Xx, zodanig dat 
y(x) 2 Yp y(x) = Y4 


waarbij Yp en Y4 voorgeschreven waarden zijn. Een dergelijk probleem 
wordt een randwaardeprobleem genoemd. Randwaardeproblemen 
komen in toepassingen zeer veelvoudig voor. 


Als uitgangspunt van de analyse nemen wij de algemene oplossing 
y(x) = A ò(x) + B blx) 


Blijkbaar moet voor een oplossing van het randwaardeprobleem 


gelden 


A ò Cx) + B vx) Yo 


A wl) + B vx) Yi 


Dit stelsel van twee lineaire algebraische vergelijkingen voor 
A en B is, in het geval dat Yp en Y4 niet beide nul zijn, een- 
duidig oplosbaar dan en slechts dan, als 

P(x) vx) 


DD = # 0 
bx) vx) 


Het blijkt aldus dat existentie en eenduidigheid van oplossingen 
van het randwaardeprobleem gebonden is aan een conditie. Om 
deze nader te interpreteren beschouwen wij nog nader een speciaal 


randwaardeprobleem, waarin wij eisen 


y(x) = 0 5 ylxy) = 0 


Men vindt nu gemakkelijk dat voor dit speciale probleem niet 
triviale oplossingen bestaan dan en slechts dan, als D = Ó. 


Wij vatten onze resultaten als volgt samen 


Definitie: Het randwaardeprobleem heet homogeen als 


y(x) = y(x) = 0. 


Stelling (Alternatief van Fredholm). 

Of het inhomogene randwaardeprobleem bezit een eenduidige op- 
lossing, of het homogene randwaardeprobleem bezit een triviale 
oplossing. 

Als voorbeeld, en ter nadere toelichting, beschouwen wij de 
eenvoudige differentiaalvergelijking 


a° | 

EL 4 Ay = 0 
2 

dx 


Ss OM>0. 


Zij OS xsS1 53 y(0) = yo» yi) = y, 


Met de eerder beschreven methode voor lineaire stelsels met 
constante co@ficienten (of eenvoudig door substitutie) vinden 


Wij dat functies 


oi) x WA 


| 


voldoen aan de differentiaalvergelijking. Doór lineair te com- 


bineren schrijven wij 


y(x) = A cos /Ax + Bsin /Ax. 
Er moet gelden 
A = Yo 


A cos /)\ + B sin /\ = Y4 





Als sin /X #0 (d.w.z. /À f mr) dan kunnen wij eenduidig op- 


lossen 
Y17Yg C9S /) 
sin /À 


B = 





Blijkbaar yy * Y4 * 0 > A = B e= 0, d.w.z. het homogene probleem 
heeft geen niet triviale oplossingen. 

Als nu /À = mr, dan is het duidelijk dat de functie B sin mix; 
met B een willekeurige constante, oplossing is van het homogene 
randwaardeprobleem. Het inhomogene randwaardeprobleem keait 
in het algemeen geen oplossingen, maar uitzonderingen kunnen 
zich voordoen, met name als toevallig Yo cos mm = Ys: 
Bijvoorbeeld, als A = m2 en Yo 5 TY4 dan voldoet de functie 
Yo COS TX 


aan alle eisen van het randwaardeprobleem. 


Echter,de functie 
Yg COS TX + B sin TX 


voldoet eveneens, en wel voor alle B. De oplossing is dus niet 
eenduidig. Men ga na dat dit in overeenstemming is met de uits 


spraken van de alternatief voor Fredholm. 





VARIATIONELE PROBLEMEN 


Inleiding. 

De benamingen honest ne tet toneis formulering", 
"variationeel principe!" etc. hebben nn Boksen in de mecha- 
nica. Men vond dat vele problemen op twee nagenoeg equivalente 
wijzen konden worden geformuleerd: door middel van differen- 
tiaalvergelijkingen of als opgave een functie te vinden die 
een bepaalde uitdrukking (een functionaal) minimaliseert. 

Voor sommige problemen was de tweede formulering zelfs natuur- 


lijker dan de eerste. Bijvoorbeeld: 


Men: verbindt twee punten A en B in het verticale vlak met een 
buis en en onder invloed van zwaartekracht EA stoffelijk 
punt van A naar B bewegen (A te hoen dn B, wrijving wordt 
verwaarloosd). Men vraagt de vorm van de buis zodanig te be= 
palen dat de beweging in de kortst mogelijke tijd plaats vindt. 
Deze opgave wordt door toepassing van de wetten der mechanica 
vertaald als volgt: 

Bepaal een differentieerbare functie x > y(x) zodanig dat 


v. 


res DE ox 
= É TIES: 
1 
een minimum aanneemt. Hierbij zijn Xx, en Xx, gegeven en wordt 


1 2 


‘geëist dat y(x) en y(x) voorgeschreven waarden aannemen. 


De moderne variatierekening houdt zich bezig met het bepalen 
van extrema van funectionalen gedefinieerd op vector ruimten. 
De ontwikkeling werd gestimuleerd door twee oorzaken: 


1. In vele problemen van de toegepaste wiskunde is de opgave 





en end ani enn eenen Tolenaars denke ad an 





om funktionalen van allerlei aard te minimalizeren de natuur- 
lijke formulering. Bijvoorbeeld in de problematiek van zoge- 
naamde optimale besturing van processen (optimal control). 

2. “Variationele aanpak!" bleek ook vruchtbaar te zijn voor pro- 


blemen geformuleerd door middel van differentiaalvergelijkingen. 


Het elementaire variatie probleem. 


In wat volgt zijn alle grootheden reêel. Zij D een gesloten inter- 
val D = {xla <x Sb}. V(D) is een nader te preciseren verzameling 


functies x *» y/@)die een lineaire ruimte vormen en die voorzien 


is van een norm. Zij nu F een gegeven functie (x,y,z) > F(x,y,2), wij 


beschouwen de uitdrukking 

Jy) = î F(x,y(x), y!'(x)) dx 
he: it element ve V(D) wordt aldus een regel getal J(y) toege- 
voegd. Wij noemen J een functionaal op V(D). 
De opgave is nu als volgt: 

Bepaal Yo E V(D) zodanig dat voor elke y € V(D) geldt 

Jy) S EAP 

De existentie van zulk een element Yo E V(D) moet uiteraard nog 
worden aangetoond. 
Eni het bovenstaande is al stilzwijgend aangenomen dat y differen- 
tieerbaar is. Voorts kunnen nog verschillende "nevenvoorwaarden'' 
op y worden opgelegd die in het nader preciseren van V tot nikes 
drukking zullen komen. Wij zullen hier zogenaamde "randvoor- 


waarden!" opleggen. 


Wij hebben al gezien dat men "randvoorwaarden" kan hebben van het 


type: y(a) en y(b) gegeven. Ook is het denkbaar dat y(a) en y'(a) 
voorgeschreven zijn, of dat y'(a) en y'(b) zijn voorgeschreven. 
In al deze gevallen kan het probleem worden gereduceerd tot 


zogenaamde "homogene!" (of 'nul') randvoorwaarden, als volgt: 





N Om de gedachten te bepalen, zij 
| | yla) = a,y(b) = 8 voorgeschreven. 
L Zij U € V een willekeurig gekozen functie zodanig dat Ula) = a 


en U(b) = 6. 


ak earteer annen 


Zij nu y(x) = y(x) - Ux). 


| D le) 
ln J(y) = Sf Flx‚y +U,y!+U!)dx = f f(x‚y,y')dx = J(y) 
da | | 


{u 





Wij zoeken Yo E V zodanig dat Yo (a) Yo (D) = 0 en 


| TG) STG) 





Û voor alle y S V waarvoor geldt y(a) vb) = 0. 
Wij kunnen dus het probleem beperken tot 
vEy Cy 


VD) = {yly E VCD); Vla) = Y(b) = 0} 





— In wat nu volgt wordt het probleem verondersteld homogeen te 





Zijn gemaakt op boven beschreven wijze, en de "“bovenstreping'" 


van alle symbolen verder achterwege gelaten. 


De eerste variatie. 





9 Wij vatten nogmaals het probleem samen: 


Bepaal Yo E V zodanig dat voor elk y € V geldt: 


Jy) S J(y) 


D 
diy) = f ECR) Vind, 
a 


De vraag is nu: hoe kan men zulk een element Yo E v opsporen? 
Wij putten de inspiratie uit een aanzienlijk eenvoudiger pro- 
bleem: het bepalen van extreme waarden van een funetie E(X) 

gedefinieerd op zeg a $x Sb. Als £ differentieerbaar is dan 


— geldt in een eenvoudig extremum £/ Rn 0. 


(Xx 




















nr rn sen Deden aande ann an madeleines eneen eneen ee rra ia eeh ae nn dn Re erheen pen. ES . 


sE 


Het opleggen van deze noodzakelijke voorwaarde leidt meestal 
tot bepalen langs eenvoudige weg van de waarden van x waar- 


voor extreme waarden kunnen optreden. 


Wij keren nu terug tot onderzoek van J(y). 


Schrijf y = Yo + mn waarbij n een willekeurig element van V 


is en T een reëel getal, zeg -T, ST S T,. Voor elke n wordt 


0 ê 


nu een reêle functie I(T) gedefinieerd door 
I(T) = J(y, + Tn) 
Kennelijk moet I(tT) een minimum bezitten voor Tt = 0. Als I(t) 


differentieerbaar is, dan geldt: 
dIi(Tt) 
S 5 nd no PA 7 5 
(Cyd SJ(y), Vy E Vh» (Sg * OsFn S Vh 


We veronderstellen nu 
F(x,y,z) is differentieerbaar naar y en z (de afge- 
geleiden duiden wij aan met Ea resp. F_). 

Wij hebben: 


Db 
JCygtn) = | F(xsygoodtin(x), Yo tin! dd dx 
a 


Voor elke n € Vn) vinden wij: 





dI(T) Pip ar Leden vite e me) 4 
dT 7 e y >Yg n 9 Yo n n 
+ P_OxygOodtin(x), ype tin! ._n!'(x)}dx 
Dus: 
ge A de HS 
dT T=0 0? 


waarbij wij als afkorting hebben ingevoerd 
_b 
! and t e 
<J (yon? Aaa zie hide yo) nx) + 


+ Fy Gs Yo (2) ° n'(x)}dx 


Men schrijft ook vaak formeel 
SJ = <d'(y den > rf 


en noemt $J "de eerste variatie van J'". 


Samenvattend hebben wij gevonden: 


Als F(x,y;z) differentieerbaar is naar y en z 


dan geldt: 
Uy) Ss J(y),vy E Vv} > 


md {<d! (yo) nes 0 Am S Wi. 
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De vergelijking van Euler. 


Door invoeren van verdere hypothesen kan het variatieprobleem 
worden teruggebracht tot een randwaardeprobleem voor een 
differentiaalvergelijking. 
Veronderstel dat 

Te Ey): Y) 
bestaat en continu is op [a;bl. Wij kunnen in êJ partieel inte- 
greren en er volgt: 

b 
<J' (yo) n> = d Ly (2) . nlx)dx, Yn E V 


waarbij: 


d 
LY 0 p Enyo) 7 dx ECXsyYpsYg) 


We veronderstellen verder, dat V de continue funktie op [a,b] bevat. 
Men maakt nu gebruik van: 


Lemma (Lagrange) 


Zij u(x) een continue functie op [a,b] en zodanig dat voor elke 


continue functie v(x) met vla) = v(b) = 0 geldt 
b | 
Sf ulx) vx) dx = 0 
a 


dan geldt: u(x) = 0, Yx € [a,bl 
Dus , Yo) is een oplossing van de differentiaalvergelijking 
Ly = 0 
(met randvoorwaarden Yola) = Yb) = 0.) 
Lyo = 0 wordt de (bij het variatieprobleem behorende) vergelijking 
van Euler genoemd. 
Wij onderzoeken nu nader de gemaakte veronderstellingen. Formeel 
differentiëren levert: 


d : 
a Eey Yo) = 











waarbij men telkens voor F moet lezen Fy, 00): VED, 


Dus, de veronderstelling dat En „ys Yp) bestaat 

en continu is op [a,b] impliceert dat F(x,y,z) aan een aantal 
differentiabiliteitseisen voldoet en dat de oplossing Yo EotT 
de klasse van twee keer continu differentieerbare functies 
behoort. 

Wij besluiten met een eenvoudig voorbeeld van de vergelijking 
van Euler. 

Lat FlRsysg) = px) z + ay + 2f(x)y waarbij ps; q en f 
gegeven continue functies zijn en p'(x) bestaat. Men vindt na 


uitwerking 


Ly, za (plx)yo) = ax) Yp „Tx. 


De methode van Ritz. 





Veelal zal het expliciet oplossen van het randwaarde probleem 
voor de Euler vergelijking buiten de mogelijkheden liggen en 
zal men zijn toevlucht moeten nemen tot benaderingsmethoden of 
numerieke procedures. 
Ook kan het zijn dat het oorspronkelijke variatie probleem geen 
oplossing bezit binnen de klasse van twee keer continu differen- 
tieerbare functies, dat F niet aan alle gestelde eisen voldoet 
en de vergelijking van Euler dus niet kunnen worden afgeleid. 
Wij keren daarom terug naar de meer fundamentele formulering 
van ons probleem d.w.z. 

J(yp) S J(y) Vy € V 
Volgens de methode van Ritz wordt nu Yo) benaderd door een 
eindige lineaire combinatie van "geschikt gekozen!" functies 
bx). Wij zullen straks ingaan op de vraag, in welke zin een 


benadering voor y‚, wordt verkregen, en ook op de vraag hoe men 
0 


de functies dx) “geschikt! moet kiezen. Eerst werken wij for- 


meel de procedure volgens Ritz uit. 


In 

Wij definieren: Re > ed) waarbij Ce, voorlopig nog 
n=1 
willekeurige constanten zijn. 


b (m) 
Beschouw: Ty) = f Fly De), W) dx 
a 


Wij wensen nu Cs n=1,...m, zodanig te bepalen dat voor de 
Cm) 


(m) 


m 


corresponderende y die wij aanduiden met Yp > geldt: 
my) ey 
Nu is Ty) = Ile,»---C…) OP te vatten als een reele functie 


van m variabelen c ce, en een noodzakelijke voorwaarde voor 


ze 


het minimum is 


Ei Erten 5 me EE 


dl 


Je: 
n 


z= 0, n = 1s.«.;m 


Aldus wordt verkregen een stelsel voor m algebraische verge- 


en} oplossing van dit stelsel bepaalt 


Be | 
Uiteraard hopen wij dat De een benadering van Yo zal 


lijkingen voor c 


(m) 
Yo n 
ZL Te 
Het is gemakkelijk na te gaan dat als ge en Be bestaan, men 


vindt 


(Cm) 


<d' (yo 


}s Ò > SN MS Anak 


Als bovendien aan de eisen gebruikt bij de afleiding van de 


vergelijking van Euler is voldaan, dan kan men werken met: 


erde 
S Ly (x) bede eN 4 me Asanes 
a 


(m) 


Nu dient nog te worden bewezen dat Yo > Yo als m > eo, Voor 
enige belangrijke klassen van problemen is dit gelukt, enkele 
resultaten zullen wij in de volgende paragraaf samenvatten. 

In toepassingen gebruikt men vaak de bovenbeschreven methode van 


benadering van Yo door ge 


zelfs in gevallen waar formele 
convergentiebewijs vooralsnog sntpteekt. Veelal is gebleken 
dat de verkregen "benadering" zelfs bij kleine waarden van m, 
de door experiment waarneembare fenomenen zeer goed weergeeft. 


(m) 


Wij gaan nu in op de vraag: in welke zin zou y, ‘ als benade- 


ring van Yo kunnen worden opgevat (d.w.z.: wat is de fout? Of, 
nog anders gezegd: als jn naar yo convergeert voor m > eo, In 


welke zin vindt de convergentie plaats?). 
Om zich een indruk te vormen over wat men kan verwachten, 
beschouw het eenvoudige probleem waarin 


b 
JY) =S Ay! GOS - 2) y 0) }ax 
a | 


Ee 


Het is eenvoudig na te gaan dat voor Yo moet gelden 


D 
S (yo) ._n!'(x) - f(x)nlx)}dx = 0 Yn E V 
a 


Volgens de methode van Ritz wordt Yo benaderd door Ere 


waarvoor geldt (m) 
b dy (x) 
SU 


a 


2 (m) 


0 (x)}dx = OQ 


- f(x)y 


Uit deze twee relaties kunnen wij afleiden dat 
b dy Cx) aye oo 2 b CD 
di en Cn F'dx S= Â fy) " Yo (x)]dx 


0 en 


Stel nu dat men zou kunnen aantonen dat Yo Yo als m » ee 


in de zin van uniforme convergentie, d.w.z. 
Cm) 
0 


Sup |y 


(A) — Va O als m > oe, 
ab 


Dan zou volgen 


b dy) aya”) 0) 


On Cek 2 > 
d nt EE }“ dx > 0 als m > oo 
ay () dy 
Men zegt dat EEn ad: van als m > ee in de zin van "kwadratisch 


gemiddelde". Dat deze convergentie niet impliceert uniforme 


convergentie, blijkt uit het nu volgende voorbeeld 





s an (m) pe 4 —_ MX < < 
41 Yg£) Yo (Cx) = ie „xl 
Wij vinden 
Sup yo - ad > 0 als m > oe 
OSxS1 
en ook (Cm) 
1 dy, (x) dy, (Xx) 9 1 he 
Á Eer NE tE dd = nn (1 -e ) > 0 als m >» oo 
Echter ES 
dy dy 
EUD le Sd 
OSS1 


Bovenstaande overwegingen maken het wellicht niet al te verras- 
send dat men, voor de verdere ontwikkeling van de theorie voor 


de ruimte V kiest een Hilbert ruimte waarin inwendig product 








gedefinieerd is door 


b 
(u,v) = f Eee eee + u(x) v(x)}dx 
a 


De norm is dan 


b 
sit = f (Eee je + LyGol Jax 
d 


Een dergelijke Hilbert ruimte wordt een Sobolev ruimte genoemd 
en aangeduid met H*(D). Daar wij geëist hebben y(a) = y(b) = 0, 
werken wij in een deelruimte Hj(D) C si 

Van een Hilbert ruimte is bekend dat er een basis in EED dane 
een aftelbare rij lineair onafhankelijke functies met de eteensenas 
dat elk element uit V met een "willekeurige nauwkeurigheid!" kan 
worden benaderd door een lineaire combinatie van de functies die 

de basis vormen (de "willekeurige nauwkeurigheid! moet uiteraard 
worden verstaan in de zin van de ingevoerde norm). 

Het ligt nu voor de hand om als "geschikte!" keuze van de functies 
bx) in de Ritz-methode te beschouwen de keuze van &'s die in 


Hoi) een basis vormen. 











Lineaire problemen. 


Een variatie probleem wordt lineair genoemd als FOXY gsYn en 





F_(x5y0»Yg) lineaire functies zijn van Yg en Yo: In dat geval 
hebben wij als algemene vorm: 


Db 
ÖJ = Öffalxdyopladm! Cx) Bd yn dB dy admi) 
da 
+yGo)ypGon bats f Cn! Ox) tga) xda 


waarbij a, B, Y, f en g gegeven functies zijn. 


Wij definieren nu een bilineaire vorm 


Lyn) z } (apr! 0018,009p GONG E0Y GOGO 
tyG)yg Gn) }dx 
en een lineaire vorm: 
b 
S(n) = =S {É£(x) n'(x) + glx) n(x)}dx 
a 
Ons variatie probleem is aldus teruggebracht tot het bepalen van 
Yo E V zodanig dat 
Lyon) z S(n) Yn € V 
De nu volgende twee resultaten (waarvan het bewijs hier achterwege 


wordt gelaten) geven de fundering van de lineaire theorie. 


Lemma (van Lax en Milgram). 


Zij V een Hilbert ruimte, Ly gn) een continue bilineaire vorm 
op V met de eigenschap dat 
Lln;n) > ollnll, o > 0 
Zij voorts S(n) een continue lineaire vorm op V. 
Dan bestaat er Één en slechts één element Yo E V waarvoor geldt 


Lyn) = S(n) Yn € V.. 


Stelling (convergentie van de Ritz-methode). 


Laat L(y,n) en S(n) voldoen aan de eisen van voorgaande’ Lemma. 





(Cm) 


Zij {é_} een basis in V en Yo 


functies gedefinieerd door 


(m) en 
Yo (xj = B 


terwijl de coëfficiënten on voldoen aan 
LYC 6) =S), n= 1 
Yo ’ De a Pe ‚DN 5 Agee sM 


dan geldt 


(Cm) 


Yo > Yo als m > ee in V. 


Toepassing. 


Beschouw het lineaire variatie probleem. 
E 2 2 
Lln,n) = f {aG)ln'Cx)l + Blxdn' Ge)ntadtyoadl nx) Jax 
a | 
Stel B(x) differentieerbaar. Door partieel integreren volgt 
5 2 2 
L(n.n) = Sf faGoln'GO + [yCx) - HB'GOEN(x)] Jax 
a 
Als nu a(x) > O0 en y(x) - 38'(x) > O0 op [a,b] dan bestaat er een 


constante T > 0 zodanig dat 


dE 


Ln.) >olnlf v ne Hij 


Het lemma en de stelling zijn nu van toepassing. 
Men overtuigt zich gemakkelijk dat de in de stelling geformuleerde 
vergelijkingen voor Co nr 1,...,m dezelfde zijn als de eerder 


afgeleide relaties 


( 


<J' (yo 


On bj == Ones Ayse m. 


Wij vinden aldus dat de Ritzbenadering convergeert in HJCD). 


Opmerking: De boven gegeven lemma en stelling gelden ook voor 


functies van meerdere variabelen. Men beschouwt dan x > Y(x), 








Pad 
El 
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xEDGC R* en de functionalen worden door herhaalde integralen 


over D gedefinieerd. 


De methode van Galerkin voor differentiaalvergelijkingen. 





Laat Ly, = 0 een differentiaalvergelijking van de tweede orde 


symboliseren, bijvoorbeeld 


Lyg =S (POOyY)) - ay, - FG) = 0 


Men wenst op interval [a,b) de functie Yg zodanig te bepalen 
dat yla) = y(b) = 0. 
Wij hebben gezien (bij de afleiding van de vergelijking van Euler) 


dat: 


Ss | 
U Lygo) node = 0, Vnev}= 
ad 


> {Lyg = 0, ypla) = y,lb) = 0} 


Wij kunnen dus, in plaats van het randwaarde probleem voor de 
differentiaalvergelijking, een variationeel probleem oplossen. 
De methode van Galerkin voor de differentiaalvergelijkingen 
bestaat uit het oplossen volgens Ritz van het corresponderend 
variationele probleem. 


Men spreekt daarom ook vaak van de methode van Ritz-Galerkin. 


Samenvattend bestaat de Galerkin methode uit het volgende: 
De oplossing Yg) van Ly = 0 met Yg la) = Yb) = 0 wordt be- 


naderd door de functie 


(m) En 
Yo (Xx) = EA et (x) 


waarbij {ò_} een "geschikt gekozen! basis is in V en een = 1,..m 


voldoen aan 





e (m) 
S Ly (xXx) « b_(x)da Es Mer beke g 
a 


Ten aanzien van de convergentie geldt dat indien p en q zodanig 


zijn dat voorgaande theorie kan worden toegepast, dan 


(1) 


(m) 
Yo ee C 


> Yo als m >» ee in H 


Wat deze convergentie kan inhouden lichten wij nogmaals toe aan de 


hand van een eerder behandeld voorbeeld. 


stel dat Yo) ie 1 Go = 2 Nkk ‚,0SxsS$<1. 


Cm) 


Wij hebben ly, Je + 0 als m > in HD). Voorts vinden wij 


Sup IyoGO - yo > 0 als me, 


O1 
(m) 


Dus: Yo mag wel in de gebruikelijke zin als een "zeer goede" 
benadering van Yo worden opgevat. 
Ten aanzien van de afgeleiden geldt 


dy, (x) ay) 





Bp I= ds 
OK dx ce MEE 
echter am) 
dy (x) dy ER) 
Sup ane an 5 ae anne > QU als m ” es 
S1 | 


voor elke ê > 0, hoe klein ook. 


(m) 





dy, | 10 
Dus: en is ook een "zeer goede! benadering van me 3 echter met 
uitzondering van een dittekourte kleine omgeving van Xx = 0, 


Of een dergelijk "locaal falen! van de benadering kan worden ge- 
accepteerd hangt af van de aard van het behandelde probleem. 

Tot slot werken wij als voorbeeld de Galerkin-methode nader uit 
voor het probleem 


ds 


f(x);s O0 SXxS1 


0 
5 aCx)yg 





dx 


Ĳ 
en, 


Ygt0) = Yo) 


Wij kiezen als basis 


Ò_ = sin nTX 
n 


en benaderen Yo door 


(m) D & : 
Yo ( = c. sin kmx 


k=0 KS 


De coefficienten Ce moeten voldoen aan 


D dee ans D 
S maen nk — a(x)y (x)} sin nmx dx TT flxjsin nnx dx 
a dx a | 


Twee keer partieel integreren levert 


b 
S n° n° - g(x)} AE sin nmx dx = = f f(x) sin nmx dx 
a a 


(m) 


0 volgt 


Door substitutie van y 


b 2 2 b 
Ek 
k=1 a a 
NEN 
Dus , Cy k = 1,...‚m worden verkregen als oplossing van een stelsel 


van m lineaire algebraïsche vergelijkingen. 


Variatie problemen en optimale besturing. 


Zoals al eerder opgemerkt, liggen vele moderne toepassingen voor 
variatierekening op het gebied van zogenaamde optimale besturing. 
Wij zullen thans een typische probleemstelling uit dit gebied van 


de toegepaste wiskunde beschrijven. 


Stel dat de toestand van een systeem voor elke x € [a,b] gekarak- 
teriseerd wordt door n functies y‚(x),...sy_ (x). 

2 n 
De funetionering van het systeem wordt beschreven door een stelsel 


differentiaalvergelijkingen 





S (n°mT°-qlx))sin knx sin nmx dx = - f f(x) sin nmx dx — 





dy; 


8; (xs; Y4>-*:Yo> U ‚>U) hk sheet 


zn 12 


waarbij u ‚U, OP [a,b] gedefinieerde en voorts nog nader te 


B 
specificeren functies zijn. 
In vector notatie schrijven wij 


ay 


a Sn 


Y worden de toestands variabelen genoemd en U de stuurvariabelen. 
Door verschillende keuzen van U te maken kunnen wij het systeem 
beïnvloeden en streven naar een "optimale!" functionering volgens 
een nader te specificren criterium. Bijvoorbeeld: Y(a) en Y(b) zijn 


voorgeschreven, U dient zodanig te worden bepaald dat de functionaal 


D 
JY Ui == f frl, Ulxdh dk 
a 
minimaal is. 
In het allereenvoudigste geval waar in n = r = 1, G(x;y;su) = u; 


vinden wij hier terug het door ons eerder beschouwde elementaire 
variatie probleem. 

In het geval n = r, G(x,Y,U) = U vinden wij een generalisatie van 
het elementaire variatie probleem tot vector-functies. 

In het algemene geval worden wij geconfronteerd met variatie pro- 


blemen van een nieuw type waaraan veel modern onderzoek is gewijd. 


Literatuur (een leidraad). 

Er bestaat een zeer omvangrijke hieu over de zogenaamde klas- 
sieke variatierekening d.w.z. variatierekening geïnspireerd door 
problemen van de fysika. Een goede inleiding vindt men in de bijdrage 
van Bottema in "Handboek der Wiskunde! deel 2, en voorts in hoofd- 
stukken b& en 6 van Courant-Hilbert "Methoden der mathematischen 


Physik I'. Men bedenke wel dat dit laatste boek, alhoewel zeer be- 





ind en nog steeds zeer veel nn geschreven is voor Ge 
moderne analyse haar huidige vorm heeft gevonden. (Zo spreekt 
Courant over "Funktionenfunktionen! wanneer hij funktionalen 

wil aanduiden). 

Een modernere inleiding vindt men in de bijdrage van Moisseev en 
Tikhomirov in "Mathematics applied to physics!" (Editor E. Rou- 
bine) Springer Verlag (1970). Deze auteurs voren niet in de Sobo- 
lev ruimte Hj(D), maar gebruiken de supremum norm. Zij behande- 
len ook niet de Ritz-Galerkin methode. Formulering en resultaten 
in Hilbert ruimte vindt men in bijdrage van Lions, eveneens in 
“Mathematics Applied to physics''. 

Uitgewerkte voorbeelden van toepassing van Ritz-Galerkin enen 
vindt men (nogal langdradig behandeld) in "“Approximate Methods 

of Higher Analysis! hoofdstuk 5 van Kantorovich en Krilov, Noord- 
hoff 1964. 

Alle bovengenoemde werken bevatten uiteraard uitvoerige verdere 


literatuur verwijzingen. 





AoXMPTOTLEK 


Ter Inleiding. 


Zij Ô(t) een reële functie gedefinieerd voor alle t € T C R , en 





zij t, een limietpunt van T. Stel dat men de functie ® wil ana- 


lyzeren in een omgeving T, van to Ten behoeve van deze analyse 
wil men &(t) vervangen door een geschikt gekozen benadering $(t). 
Een dergelijke benadering wordt "asymptotisch' genoemd indien de 
benadering beter wordt (in een nader te preciseren zin) naarmate 


t dichter bij t, ligt. 


0 
Als een eenvoudig voorbeeld beschouw 


Ë 


DCE) = Ct) + e ò,Ct) 


waarbij dt) en Ò,‚Ct) rationale functies zijn, gedefinieerd op 
O & t < eo, terwijl bct) begrensd is op dit interval. Stel dat 
men geinteresseerd is in d(t) voor zeer grote waarden van t. 


Het is duidelijk dat 


lect) - 6,CH] Sc e* 


\ 


waarbij c een constante is. Dus, voor grote waarden van t, kan 


de functie bt) als asymptotische benadering van ®(t) dienen. 


Als tweede, en sterker sprekend voorbeeld, beschouw 


Z 


Dt) = 25 


® 


eeN 


Ss tE 2 


om 8 
et 
N 
ev 


Men kan aantonen dat 


oo k 
Ö(t) = e {-Int - yy + TZ (-1)K71 Ln 
k=1 k k; 








nnn ee ee etn kn ee. rr ennen nak 
Ge eene 


De machtreeks in het rechterlid convergeert voor alle waarden 
van t > 0. Desondanks is deze formule voor het onderzoek van 
b(t) voor grote waarden van t volkomen ongeschikt. Wil men nume- 
rieke waarde van ®(t) voor een grote waarde van t met een rede- 
lijke nauwkeurigheid bepalen, dan is een zeer groot aantal 
termen van de machtreeks nodig. 

Bij een voorgeschreven nauwkeurigheid neemt het aantal benodigde 
termen toe naarmate t toeneemt. 

Een eenvoudige asymptotische benadering van Ò(t) voor grote 


waarden van t kan als volgt worden verkregen: 


Ae En ke E _ 2 
Schrijf: ke 1 + [ST il = 1 ET 


Substitutie in de integraal die ®(t) definieert levert 





oo oo 
bt) =f eZ az -f ef Ar dz 
0 0 ai 
Zl 
zi + R 
waarbij 
— 1 
IR‚l =S 82 Ar az < f e°7 zdz = So 
0 0 t 
Dus: d(t) = Ll kan dienen als asymptotische benadering van Ô(t). 


Ë 
De nauwkeurigheid wordt geven door de schatting 
1 
oct) - HH] < 5 
t 


Men kan de asymptotische benadering verder "verbeteren!" door 


als volgt te werk te gaan: 


pe 





m m 
Schrijf: a * E (-z) SE (A Bla 
k=1 “k=1 
m m 
B k-1 (C-z) 
= 5 (z) A 


vl 
eN 


Substitutie in de integraal levert nu 





d(t) = E Ci Eke)! t7K « Ro 
k=1 
waarbij 
oo eo 
|R | Á Ee AT dz < / et2 2Fgz = m! mri 


Wij merken op dat 


R_ > 0 als t * 
m 


en zelfs 
R__ 
A > 0) als t *ee, 
m 
En k-1 -k 
Echter, de formele reeks £ (=1) (k-1)! t is divergent! 
k=1 


Men kan dus niet verwachten dat men de nauwkeurigheid van een 
asymptotische benadering, bij een vaste t (hoe groot ook) wille- 
keurig kan opvoeren door meer termen van de reeks te nemen. 

Bij het hanteren van asymptotische benaderingen is een flinke 


mate van begrip en voorzichtigheid geboden. 


Asymptotische grondbegrippen. 








Zij nu ®(t) en Y(t) twee functies, gedefinieerd voor t € T C R 


a er Ee hi na EEE ren TE ve lieke Srei hehert asch hr Be ndnanlnhein” zhe 


en z1j to een limietpunt van T. Om het gedrag van ® en Y, als 


tE tE te vergelijken gebruiken wij de volgende orde-symbolen. 


oke 
Definitie. ® = O(Y) voor t * to indien er een constante k bestaat 


en een omgeving U van t‚ zodanig dat 


0 
lol SS k|Y| voor t € U 
® = o(Y) voor t * to indien 
lim 7 z= 0. 
ttg 
Door toepassing van orde-symbolen kunnen functies asymptotisch 
worden geordend. In alles wat nu volgt bestuderen wij (zonder 
het iedere keer opnieuw expliciet te vermelden) het gedrag voor 
t > to» waarbij ty een limiet punt is van het definitie gebied 


der te beschouwen functies. 


Definitie. Een rij (6, Ct} k = 1, 2,... is een asymptotisch ge- 
ordende rij, of kortweg een asymptotische rij, indien 
voor alle k geldt 


Ori 5 OlD) 
m 


Een reeks ZE adt), waarbij a, (numerieke) constanten 


k 
k=1 
zijn, is een asymptotische reeks, indien Kn, een asymp- 


totische rij is. 


Met behulp van asymptotische reeksen kunnen wij op een systematische 


manier "alsmaar betere" asymptotische benaderingen definiëren. 


Definitie. Zij gegeven een functie ®(t). 6,Ct) is een asymptotische 


benadering van ®(t) indien geldt 


Ò -— Ô, = o(s,)- 








| m 
De asymptotische reeks 2 a, Pr 
k=1 B 


ontwikkeling van &(t) indien geldt 


(t) is een asymptotische 


a 
1 


@ 
| 
u M3 


ië k?k à o(Ô) 


Opmerkingen. 


m | NE: 
Als 2 adt) een asymptotische ontwikkeling van ®(t) is, dan 
k=1 Ek 

schrijft men vaak 

m 

blt) — Za dt, (t) 

isd kk 
Als een asymptotische ontwikkeling voor alle m geldt dan schrijft 
men ook 


b(t) > Za dt) 
k=1 kk 


Deze schrijfwijze is symbolisch, daar de reeks ook divergent kan zijn. 


Wil men aangeven dat ®(t) begrensd is, dan is het gebruikelijk te 
schrijven 
b(t) = O(1) 


Men schrijft ook &(t) = o(1) om aan te geven dat Ò(t) * 0 als t > to: 


Rekenen met orde symbolen. 


De nu volgende relaties kunnen zeer eenvoudig worden bewezen: 
1. Old) - OW) Od.) 
old.) 


3. Als & = Olp) dan is: 


2. Old) - olk) 


3.1. Od) + Olp) = OCW) 
3.2. o(&) + olp) = ol) 
3.3. Od) + ol) = OCW) 





Het hanteren van asymptotische ontwikkelingen. 


Bij het analyseren van functies door middel van asymptotische 
benen dient men de vraag te stellen: hoeveel termen van 
de reeks moet men "meenemen" ei wat is de verkregen nauwkeurigheid. 
Bij het beantwoorden van deze vraag zijn de nu volgende overwegin- 


gen van nut. 


I. Een asymptotische reeks kan convergent of divergent zijn. Is de 
reeks convergent dan behoeft de som van de reeks niet gelijk te 


zijn aan de ontwikkelde functie. Als voorbeeld beschouw 


a 


EEn 
b(t) = Tet + e 


Voor t * ee hebben wij 


oo 
blt) — E (-1)k dr #00 
k=0 t 


De benadering door een eindig aantal termen van de reeks wordt 
dus alsmaar beter naarmate men meer termen meeneemt. 
Echter de fout kan niet willekeurig klein worden gemaakt door 


voldoende groot aantal termen te nemen. 


II. Bij een divergente asymptotische reeks kan een te groot aantal 
termen onjuiste resultaten geven. Dit blijkt uit het nu volgende 
resultaat: | 


Stelling. Zij ®(t) een begrensde functie en voor t > to 


(t) + R‚Ct) 


m 
Ò(t) = 5 a 


k=1 
met 


Orar 5 CO25 Rm sol) 


Voorts [R_(t)| * ce als m *e Dan geldt: 


Ek. 


fa 











Bewijs 


1. Voor elke vaste waarde t = t4 bestaat er een natuur- 
lijk getal u zodanig dat voor alle m # u 
& ® 
[R Ct) [R_Ct4)| 


2. Als t‚ * to dan u * ee (en uiteraard ook R(t,) * 0). 


is nagenoeg triviaal. Het bestaan van een waarde m = u 
zoals in de stelling bedoeld volgt uit |R_(t)| + ee voor 


m > eo, Het feit dat u * ee als t, zá to volgt uit 


K 
Ei z o(1). 
u 





Asymptotiek van integralen 





In toepassingen is men vaak geconfronteerd met de opgave een 


integraal van het volgende type te berekenen: 


tw(z) 


HE) = e TLZ 


Rs DW 


Voorlopig worden alle grootheden reëel verondersteld, t is een 
parameter die grote waarden aanneemt. Een eenvoudig voorbeeld 
van een dergelijke integraal (met zijn asymptotische ontwikke- 


ling) hebben wij in de inleiding reeds gezien. Beschouw opnieuw 


„Gets de 
GES WE ZT 


Wij hebben aangetoond dat 


kT k 


m Kl 
en as | 
D(tT) — kei! 1 (k-1)!t 


De behandeling in de inleiding was echter geheel zonder enige 
motivering. Thans willen wij een gedachtengang ontwikkelen die 
ook in het meer algemene geval ons in staat zal stellen asymptor- 
tische ontwikkeling van integralen te bepalen. 


Beschouw: 
oo _tz 
BLE) = | 8 f(z)dz 
0 


Dergelijke integralen worden Laplace integralen benee Als t 
groot is, dan is e EZ voor z >0 zeer klein. Men kan dus 
“verwachten! dat de grootste bijdrage tot de integraal geleverd 
wordt door nn omgeving" van z = 0. In die omgeving kan men 
wellicht f(z) geschikt benaderen en aldus de asymptotische ont- 


wikkeling van Ö(t) bepalen. 





Bovenstaande redenering leidt tot het volgende resultaat: 


Stelling: Zij 


Bewijs: 


Bj fe tZe(maz 


waarbij f(z) een begrensde functie is op 0 Sz <eoe, 


terwijl de reeks 


pag n 
fuh == EE az 
h=0 n 


convergent is voor 0Sz Sz ‚ met zo >0. Dan geldt 


voor too en willekeurige m: 


m 
ö a n 


mn 
b(t) > LL a 
h=0 o nn 


oo am 
f zie tZgz = 
No n 
Wij schrijven 

6 tz 0 tz 
Ò(t) = { e f(z)dz + fe f(z)dz 

6 

waarbij TC een willekeurige maar vast gekozen getal is, 
zodanig dat 
Oc So 
Wij weten dat er een constante M bestaat zodanig dat 


[£Cz)| <M 


Er volgt dat 
rt M _=Et -1 =-Gt 
[fe f(z)dzl < e = Olt "e ) 
Ee a 
Dus: 


ap ee 
Ò(t) = Á etZf(zjdz + oct ter tt) 


Schrijf nu 


in n 
fz) o 5 az = ff ({z) 
n=0 n m 


Uit de eigenschappen van convergente machtreeksen volgt 
dat er een constante M' bestaat zodanig dat, voor 0OSz<z, 


|£_C2)| < M'zm*! 


Wij schrijven 


m m 5 ie A a 
DAE) Se EE A zie Hodr +ff (ze Ke + Ot le t) 
n= n o Mm 


| os IN 


Expliciete berekening van de integralen levert: 


G n n-1 n=2 
f zhe tZaz = -[e tZ{Z- + DZ + aln-1)z î 
0 E +? +3 
G 
nz n: __n! 1 —=Et 
+ + En + nt) | = ‚+1 + ri e g(5st) 


0 


waarbij 8, (Et) uniform begrensd is voor OS{S$zo en t > to > 0. 


Op analoge wijze vinden wij 


m+ 2 


6 n 6 5 | N 
LE {ede taz < MJ En Clas £ Mrt 10: + 5 Ct, 
De 0 + t 


Er volgt aldus: 


' 


m 1 -Gt) 


D(t) = aad dte Wk tt e 
n=0 n ‚nti 


Het is eenvoudig in te zien dat indien er voor een functie 


v(t) geldt 


1 "Gt, 


biel = Ott “e T>0 


m= 2 


er ook geldt W(t) = o(tT FTS) en zelfs (tj = olt B 


met m willekeurig. 





ve nea nnen and kt 
Ee 5 Sn Ee 





Hiermee vinden wij 


m ' 
_ aj | n= 2} 


hetgeen de stelling bewijst. 


Opmerkingen: 

De ontwikkeling van &(t) voor het gevat dat f(z) = GED 
welke in de inleiding was gegeven, is in feite een eenvoudige 
toepassing van de nu bewezen stelling. De stelling kan nog 
worden gegeneraliseerd tot het geval dat in 0O$Sz$<zo;, Zo 0; 


geldt 


fíz) = P(z) ZE, az 


waarbij P(z) een gegeven functie is. De asymptotische ontwik- 


keling is daar gegeven door 


m ee =t 
OLE EK ad P(z)zie Zaz 
n=0 Ao 


Deze generalisatie heeft belangrijke toepassingen, bijvoor- 


beeld het geval 


f(z) = 


z z+1 


Ns 


Ontwikkelingen door partiële integratie 


Het kan voorkomen dat f(z) geen convergente machtreeksontwikke- 
ling bezit in de omgeving van z = 0. Men kan dan gebruik maken 
van een andere methode (welke overigens, in het geval dat f(z) - 
wel een machtreeksontwikkeling bezit, equivalent resultaat levert 


met de reeds behandelde stelling). 





Stelling: 


Bewijs: 


„48 = 


"E48 f(z)dz 


le ©) 
Zaj Olt) = Á e 
waarbij f(z) uniform begrensd is op O$z<oo, Stel dat 


de afgeleiden 


bestaan en uniform begrensd zijn voor 0O$zS$z2zo, Zo 0, 


en n == 1323«eesmt1. 
Dan geldt 
En | 
bl rd os Po) = eco). 
ns tn+ 


Wij sehrijven opnieuw 


LZ EZ 


f(z)dz + fe 


G 
blt) fe 
4 T 


f(z)dz 


met 0 S<t<zo. Wij weten al dat 


Nad 4 


Í se eid = Olt Te ) 
E 


De eerste integraal wordt partieel m-keer geïntegreerd 


als VoLEt: 

je -tz 1 =Ct 1 je „tz (1) 

{ e lZidz = RD, —= e £(T)} + je £ (z)dz 
Ee | n 

zfe ze CD zaz = Bte ooo EE der 


+ 


En 
Efe NZ 2 zaz enz. 


+2 


De resultaten samenvattend vindt men 


m n _ n 
dt) = ZE en £ 0) OE 2) + Ot De Ge, 
n=0 ‚nti | | 


Hiermee is de stelling bewezen. 





Kritieke punten en de methode van Laplace 


In een verdere uitbreiding beschouw nu 


tw(z) 


B 
Ò(t) = fe f(z)dz 


e 


waarbij w(z) en f(z) reële functies zijn. 

EDE mogen verwachten dat de grootste bijdrage tot de integraal 

komt uit de omgeving van de punten van het interval a<z <8, 

waarvoor w(z) een maximum bereikt. Deze punten wonen kritieke 

punten genoemd. Volgens de methode van Laplace wordt de asympto- 

tische ontwikkeling van Ò(t) bepaald door beide functies w(z) | 

en f(z) in de omgeving van de kritieke punten geschikt te benade- 

ren. dere uitwerking kan men vinden in 

A. Erdélyi, “Asymptotic Expansions! , ever iose 

H.A. Lauwerier, "Asymptotic Expansions'', Mathematical Centre Traets,1966 


N.G. de Bruijn, "Asymptotic methods in Analysis", North Holland, 1958 


Integralen met snel. oscillerende integranden 


Voorgaande theorie kan worden uitgebreid tot het geval dat tw(z) 


complex is. Wij beperken ons hier tot het eenvoudige geval 


Ò(t) 


B . 
z= f ela rbde 
a 


waarbij t, f(z) en z reëel zijn. 
Wegens de snelle oscillaties van de integrant moet men verwachten 
dat bijdragen tot de integraal van opeenvolgende deelintervallen 


van de lengte 7 nagenoeg tegenelkaar zullen wegvallen. Het blijkt 


dan ook dat een significante bijdrage uitsluitend wordt geleverd 
door omgevingen van de uiteinden van het integratie interval, 


Z == a en zZz = B. 


Stelling: Zij Ot) = 


Bewijs: 


sE 


ee f(z)dz 


2 WW 


waarbij f(z) een mt+1 keer continu differentieerbare 


functie is op aSzS<B. 


n 
zij Aaf Pas ele = ete 


dz 
Dan geldt 
m k+1 : ì | | 
on — r Ee telt sl pj-eltt £ MO Cap}. 
k=l (it) 


Volgt door partiële integratie. 





Slotopmerkingen 


In het voorgaande zijn enkele elementaire grondbeginselen en 





toepassingen van de asymptotiek behandeld. Moderne ontwikkeling 


heeft vooral betrekking op functies van meerdere variabelen in 


verband met de zogenaamde storingsproblemen. 


A1 Au Eee pete variabelen en € een Rene En 
(men gebruikt e in plaats van t omdat de problemen meestal zo 
zijn geformuleerd dat de parameter kleine waarden aanneemt). 
Zij D(Xys-ee XE) een functie east intedd als oplossing van 


een differentiaalvergelijking 
L_® z= 0 voor Cajoree sk) E D 


Hierbij treedt e op als een parameter in de operator L,- Voorts 
worden voorwaarden op de rand van het definitie gebied D. opge- 


legd. 


Men wenst nu asymptotische benaderingen van ® te construeren, 
geldig voor e > 0. Het probleem is dan meestal zodanig gecompli= 


ceerd dat constructie van een expliciete representatie van Ò 


buiten de praktische mogelijkheden ligt. 


Voor de behandeling van dergelijke problemen moeten in de eerste 
plaats de grondbegrippen van de asymptotiek nader worden aangepast. 
De orde van grootte van een functie wordt gemeten door middel gas 
een geschikt gekozen norm, bijvoorbeeld 
|l = max |ò| 
D 
Men merke op dat een functie van verschillende orde van grootte 


kan zijn in verschillende deelverzamelingen van D. 








De definitie van asymptotische ontwikkelingen brengt nu ver- 
gaande moeilijkheden met zich mee. Men kan, in analogie met 


het voorgaande, invoeren de asymptotische reeksen 


m 
Ù a, (xjseeesX, db le) 8 Ö+4 = OCH) 


k=1 


en daarmee als asymptotische ontwikkeling definiëren 


| n | 
D(xgoere XE) = en a, (Xys-eesX db le) + OCH) 


Een dergelijke ontwikkeling in het geval Ò e“ wordt een 


Poincaré ontwikkeling genoemd. 


Het blijkt echter dat in vele door differentiaalvergelijkingen 
gedefinieerde problemen de functie ® geen ontwikkeling van 
boven gegeven vorm in D bezit. Dergelijke problemen worden 
singulier genoemd en zijn onderwerp van vele recente onderzoe- 
kingen. Een inleiding in de boven geschetste problematiek kan 
men vinden in: 


W. Eekhaus, "Matched asymptotic expansions and singular 


perturbations!, North-Holland, 1973. 


Praktikum Toegepaste Analyse 

| Differentiaalvergelijkingen 

1. Gegeven is de functie f(x) op (a,b). Bewijs dat als f(x) kont inu 
differentieerbaar is op (a,b), f(x) Lipschitz-kontinu is op 


[e‚d] C (a,b). 


2. Gegeven is de vergelijking 


dy 
dx 


y log y» y > 0 

= 0 ss y = 0 
We zijn geïnteresseerd in oplossingen in het x,y-vlak die gaan door 
het punt Xyp Yp In welke punten van het x;,y=vlak is voor deze 
oplossingen niet aan de voorwaarden van de existentie- en eendui- 


digheidsstelling voldaan? Geef de oplossingen met beginvoorwaarden 


y(0) = 0, y(0) = 1/2, y(0) = 1. 


3. Gegeven is de vergelijking 

dy … 2 ) e= 

mr 1 +y ; y(x) 0 
Xo E RÎ 


a. Ga na dat voor dit probleem aan de voorwaarden van de existentie 
en eenduidigheidsstelling is voldaan. 


b. Hoe groot is het interval waarop de oplossing bestaat? 


hb, In modellen voor de groei van de wereldbevolking (of de groei van 
een deel daarvan, of de bevolkingsgroei van een andere diersoort) 
geeft y(t) het aantal personen aan op tijdstip t; yo) = Yo 0) en 
het geboorteoverschot is r $% ber jaar. Dit percentage verandert in 
de loop der tijden kontinu. 
a. Bedenk (sen de laatste opmerking een toelichting. 


b. Een mathematisch model voor dit proces wordt gegeven door 


dy 
d t . 
Ga na dat het model een oplossing bezit die eenduidig is. 


= r(t)y, y(0) = Yo: 


c. Bereken deze oplossing. 








2. 


d., Leid een beteke af voor de tijd T waarin de bevolking in 

: aantal verdubbelt. Merk op dat deze verdubbelingstijd T onaf- 
hankelijk is dan de beginvoorwaarde. | 

e. In het rapport van de club van Rome (uitgave Het Spectrum, 
pag. 38) wordt opgemerkt dat in het jaar 1650 de wereldbevolking 
0.5 miljard mensen bedroeg, het geboorteoverschot was 0.3 % per 
aa: Neem aan dat dit konstant is, hoe groot wordt dan T? 

f. In 1970 telde de wereldbevolking 3.6 miljard mensen en bedroeg 
het geboorteoverschot 2.1 %. Neem aan dat het geboorteoverschot 
lineair is toegenomen van 1650 tot 1970. Bereken opnieuw de 
bevolkingsgroei y(t) en de verdubbelingstijd T. 

g., De aanname van lineaire toename van het geboorteoverschot 1s 
kennelijk te grof (waarom? ). Een realistischer model krijgen 
ge als volgt: Neem aan er is een tijdstip tT, 1650 < T < 1970, 
zodanig dat op (1650, tT] r konstant is en op [tT,‚ 1970] lineair 


toeneemt. Bepaal T; hoe groot wordt T in 1970? 


5. Een zeer eenvoudige theorie voor de beschrijving van epidemieën 
is de volgende. Webendessten dat er n leden in een afgesloten 
gemeenschap zijn, waarvan er p geïnfecteerd zijn ak q niet ge- 
Bikesctelsd. maar wel vatbaar; p+q = n. Voer in het relatieve 

| aantal zieke individuen x = È en gezonde individuen y = d, Laat 

n zo groot zijn, dat in benadering x en y als kontinue varia- 

belen mogen worden behandeld. De verspreiding van de ziekte wordt 

gegeven door 5%; de verspreiding is evenredig met het aantal 

kontakten tussen zieke en gezonde, individuen. 

a. Stel de vergelijking op die de verspreiding van de ziekte be- 
schrijft. 


b. Los de vergelijking op; bepaal lim x(t) als deze bestaat; 





hangt deze limiet van de beginvoorwaarden af? 





Ne ee 
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c. Noem een aantal redenen waarom het model niet realistisch is; 
probeer het model te verbeteren. 
d. Welk epidemisch verschijnsel korrespondeert met de vergelijking 


dx(t) _ ú n ee 
art nike Bx (t Ene x(t=-1)) 


‚ We beschouwen een meer met konstante inhoud V kem® waarin zich 


een hoeveelheid chemisch afval bevindt. Het afval is voortdurend 


gelijk verdeeld over het hele meer; op tijdstip t is de concen- 
tratie c(t) (dus 0 S ct) < 1). Veronderstel verder dat een ver- 
vuilde rivier in het meer uitkomt; de concentratie van chemisch 
afval in de rivier is konstant K en per tijdseenheid komt een 
hoeveelheid s rivierwater in het meer, veul eenzelfde hinveetis 
heid water per tijdseenheid uit het meer stroomt. 
a. Stel de differentiaalvergelijking op die het gedrag van c(t) 
met de tijd beschrijft en los deze op; bereken lim c(t). 
b. Los dit probleem op indien bovendien per ET een ge- 
wichtshoeveelheid chemisch afval G (de G van Geheim) direkt 
in het meer wordt geloosd; bereken lim c(t) indien deze bestaat. 


too | 
c.‚ Indien de tijdseenheid 1 jaar is, laat dan 7 = 15 een voor=- 
1 
zoo Sen voorbeeld van 
langzame doorstroming. We nemen verder aan dat de vervuiling 


beeld van snelle doorstroming zijn; 7 = 


van de rivier beëindigd wordt en de lozingen gestaakt op het 
tijdstip dat c(t) = 100 c(0) (dus aanname 0 < c (0) < ToT) 
Hoeveel tijd verloopt in beide gevallen voordat de concentratie 


c(t) weer de oorspronkelijke waarde c (0) bereikt. 


Opmerking: 


In dit model zijn een groot aantal faktoren verwaarloosd, zoals 
bijvoorbeeld absorptie, verdamping, lokale concentratie van ver- 


ontreiniging enz. 





be Praktikum Toegepaste Analyse Inleiding differentiaalvergelijkingen 
: 7. Bepaal de oplossingen van de volgende vergelijkingen door de 


punten esn) waarin existentie en eenduidigheid is gegarandeerd. 


_2 
f a dy YY 
| ° dx Xx+ 1 
B. xy xy = - Zy +i) 


8. Gegeven is de vergelijking 


dy - 
Ee f(y,x) 


Op het gebied D iniR? (voor het gemak beschouwen we skalaire 
grootheden, de overgang op vektorfuncties geeft geen essentiële 
moeilijkheden) 5 f(y‚x) is kontinu in x en y op D en Lipschitz- 


kontinu in y (Lipschitzkanstante A) ; verder geldt op D 


D |ÉCx‚y)| S M (konstante) 


(n) 


a. Als Je 


1 € Den y x) de n-de approximatie voorstelt 


de 40 
van de oplossing door Ea volgens de methode van Picard- 


Lindelöf, bewijs dan dat geldt voor de oplossing y(x) 


+1 
Eik 


A | x-Xg | 
e 





| y(x) -— vG SMA 
Ak | (n+1)! 


b. Beschouw als toepassing de vergelijking End = x° + y° met 
x EE [o,i]l en Eet . We willen de oplossing y(x) approximeren 
bijvoorbeeld z.d.d. 


( 





ly 60) = y en < € 
In 


Merk op dat we hiertoe D nader moeten bepalen, bijv. door 
Iyl <p (positieve konstante). Laat zien hoe voor gegeven 


e‚ n afhangt van p ; substitueer ter illustaratie eens p = 1 
id 


EEE ETE Ee BETE Tike dk ii a Sk hi a 
Ld 


en € = 10 en bereken dan n(e = 2.7183). 


CE EET Ee it ik ad TE ti iS ee a ik 

















10. 


11: 


Gegeven is de vergelijking met beginwaarde 


Z- Vh Rs PCOS G 


a. Benader de oplossing met behulp van de methode van gukseastene: 
approximaties (iteratie-methode van Pieandelfneerdf). 
Beschouw de beetmede als 1e approximatie, ga tot en met 
ge 

b. Geef de exacte oplossing en vergelijk deze met de benaderde 


door Taylorreeksontwikkeling. 


Gegeven is 
Vetex 0, y= 1 


a. Benader de oplossing met behulp van de methode van suksessieve 
approximaties 5; ga tot en met y(3 
b. Geef de exacte oplossing en vergelijk deze met de benaderde 


door Taylorreeksontwikkeling. 


Toegepaste Analyse. vervolg differentiaalvergelijkingen. 


Beschouw de vergelijking van de harmonische oscillator: 


Es 2 | 
De t a Xx == 0 met @ > Q 


dt 
x(0) = Xp 5 x Ìis de plaatscoördinaat 
dx … d Tr | 
ar 0) = Vo 3 TE z= V ls de snelheid. 
a. De Iste Picard geïtereerde wordt x (1) = Xp vei) = Vor 


| mn je 
Bepaal de n-de Picard geïtereerde zend, vn, 


De. Maak een schets in het fasevlak xev] van de banen corres- 


ponderend met: de exacte oplossing, de 2de, 3de en ude Pisa 


geïtereerden. 
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13, 


1u. 


he 


Merk op dat voor grote tijden (t > @) de banen van de Picard 
geïtereerden willekeurig ver afwijken van de exacte baan. 


Gegeven is de vergelijking ah = Alx)y 
waarbij A een kontinue n x n matrix is, gedefinieerd op het 
segment 1. Bewijs dat als voor een oplossing geldt 

FG) = 5 
met Xx, E I, dat deze oplossing identiek nul is op LI. 


Bewijs de stelling gegeven op pag. 11 van het diktaat. 


Teken het richtingsveld behorend bij de vergelijkingen, gegeven 
in vraagstuk 7. 
Kontroleer de verkregen resultaten met behulp van de onder 7 


uitgevoerde berekeningen. 


Beschouw op [0,1] de differentiaalvergelijking 


2 

SY + U - 0 
dx En 

met achtereenvolgens de randvoorwaarden 
a. y(0) = a, y(1) = B 


Be MEO md) B 


en Py(0) sk qy'(0) z A, Py(1) + sy!(1) = B 
Bestudeer de existentie- en eenduidigheid van de oplossingen van 
deze 3 randwaardeproblemen voor alle mogelijke reêle waarden van 


de konstanten a, 6, Ps da Pa Sa 


Î 
| 
ä 

i 
Î 
: 
Í 
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16. 


Re 


18. 


Variationele problemen. 


Beschouw de volgende funktionaal op C[ 0,1] 


ED 
J(y) = f y°(x)dx) 
Ò 


a. Laat zien, dat er precies één Yo E C[O,1l bestaat die J minimali- 
seert . Welke Yo is dit ? 
b. Laat zien, dat er géén funktie in C[O,1l bestaat, die J minimali- 


seert en tevens voldoet aan de randvoorwaarde y(0) = 4 # 0. 


In het diktaat (pag.24) wordt de Euler-vergelijking afgeleid 


onder de veronderstelling 
V2Ce={yEcla,bl ; yla) = y(b) = 0}. 


Deze eis is (veel) te sterk. 
Ga na, dat deze eis verzwakt kan worden tot bijvoorbeeld 


Vv CC en V ligt in de sup-norm dicht in oe 


Gegeven zij dat de gevallen V = C* [a,bl ; n > 0 hieronder vallen. 


a. Laat zien dat : 


N Ì 
2, + (y') Tax je 


Ch en 


iylh =tf Ly 


(1) 


een norm definieert op C [0,1] afkomstig van een inprodukt 
b-/We zoeken het element y € [0,1] met de kleinste norm, 


zoals in a gedefinieerd, dat tevens voldoet aan 
y(0) = a 5 y(1) = B. 


Los daartoe de Euler-vergelijking op. 
Waarom is deze oplossing uniek? 


Noem deze oplossing Yp” 


c. Bewijs, dat er géén element Yo e CD oa bestaat die aan de 





De randvoorwaarden uit b voldoet en Dy} S IÍvali. 


19. Gegeven is de funktionaal J,. die geminimaliseerd wordt op de ruimte 
Vv, 5 gegeven is verder de funktionaal J, die geminimaliseerd wordt 
op de ruimte V,. Nu geldt V, C V, en J, = J, op V‚- 


Bewering : y, € V, minimaliseert J, op V‚, > yo minimaliseert J, op Vv, 


Is deze bewering juist ? 


| ZU. Een klant bestelt op t = 0 bij een fabrikant een hoeveelhete 
goederen groot G. De goederen moeten worden geleverd op het tijdstip 
t= T (> 0). 

De fabrikant produceert deze goederen tussen t = 0 en t = T en 
streeft daarbij naar minimalisatie van kosten. 


ik y(t) de/hoeveelheid tot op tijdstip t geproduceerde goederen. 





Er zijn 2 soorten kosten 





T 
n n fl _ d 
| ls produktiekosten : K, = / k, GED dt 
| r 
Î ii. opslagkosten Ks ‚ k‚(y) dt 
| 
| 
Gegeven is 
k,(x) = x +} x° voor x > 0 
k,ly) = y voor y > 0 


a. Interpreteer de uitdrukkingen voor k en k, 








Î 
" 
| 
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Stel een bijbehorend minimaliseringsprobleem op (Funktionaal, 


randvoorwaarden, funktieruimte Vv). 


Los dit probleem op door gebruik te maken van de Euler-vergelijking 
en OPE 19. Onder welke voorwaarde voor G en T levert deze 
oplossingsmethode een korrekt antwoord ? 

Kunt u de oplossing van het probleem vinden, indien aan de taatdee 


voorwaarde niet is voldaan ? 


Een lichtstraal beweegt door een blok materiaal van 


A = Fo,1) naar B = EN, 
A 


L 


' 
JA 
B … 

ne 
O X B | 
De brekings index van het materiaal is n‚deze hangt i.h.a. af van 


de positie, dus n = n(x,y). De snelheid van het licht is e/n. 


. Het licht plant zich zodanig voort, dat de tijd benodigd om van A 


|» 


E 


naar B te komen wordt geminimaliseerd (natuurwet). 


Beschrijf de baan van de lichtstraal door : x,y) 1, xe [Ö, 1 


Laat zien, dat dit probleem zich laat formuleren als 


| 
minimaliseer  : À { n(x,y) „ M1+lyT)7 ax 
( 


randvoorwaarden : y(0) z 1 5; y(1) = 0 


Geef de oplossing als n constant is. 


n, voor Xx & 


Ne 


Geef de oplossing als n = n, voor x > 3 


me t constanten no en n, 


Maak hiervan een schets. 


n, sin(r) 2 


Ziet u de'wet van Snellius’ : ng Sin(i) s 





10. 


d. Geef de oplossing als n= a+gBx 
a1s;sa+B>i1 

8. Geef de oplossing als n= y+êy 
Yel;sy+6>1 


Hint : beschouw niet y(x) maar x(y) ! 


22. Het probleem op blz. 19 van het diktaat luidt : minimaliseer 


1 2 
S ‚a dx 


0 


y(0) = 0 5 y(1) = 1 


ol x 
nn 
on 





x=1 
a. Leid de vorm van dit minimaliseringsprobleem af 
b. Bereken de oplossing 


Hint : zie vr. 21e 


23. Beschouw een kontinu differentieerbare f(x) in IR, periodiek met 
periode 21. We willen f(x) op [0,21] zo goed mogelijk benaderen 


met een funktie van de gedaante 
y(x) = A, + Az COS X.+ A, sin Xx 
a. Kies als norm 
[lyll = [/ y?lxddx 12° 


formuleer een minimaliseringsprobleem en los dit op. 


ds 


b. Doe hetzelfde met als norm 


2 8 
ly = fy?) + y'20) ax 12 
0 





S- Wat vindt u als u f in de norm uit opg. a zo goed mogelijk wil 


benaderen met een functie van de gedaante 





YEx) = A, + 


n {A. cos(kx) + Len sin(kx)} 


s1 k 


AM 


Merk op, dat het antwoord niet van n afhangt. 


En eend ak Md a ae eer ae nan 


De gevonden A,, Aj» Ag heten Fouriercoefficiënten van f. 


2, Beschouw de vgl. + randvoorw. 
+ Ss Xx == f(t) 
x(0) =xX(N = 0 met s € R. 
a. Laat zien, dat deze vergelijking ontstaat als Euler-vergelijking 


van het funktionaal 


Ì id . 
J(x) = f { - Fx)? + 5x? - f(t). x } dt 
0 


Dit funktionaal wordt echter niet noodzakelijk geminimaliseerd 


of semaximaliseserd 5 


Geef bijv. in het geval f = 0 ens > 1? 


(1) een rij funkties z, met z,0) = z,(1) 0 en dz) * co 


voor n ” @ 


(ii) een rij funkties z, met z,(0) zi) 0 en J(z) > — @ 


voor n ”» ow, 


2 


Gegeven zij nu, dat s # k?. mT? voor k = 1,2,3,.... 


Voer het procédé van Ritz-Galerkin uit met als basis En 
(n) _ 


[7 


sin (k T t 


D.w.z. stel inde n=de iteratieslag X 


= M5 


a, &, 
ze Ek 
en bepaal de optimale set { hl ks dsvaden n } voor J(x 


Merk op dat in dit speciale geval an) niet van n afhangt. 


mo 


[2 


Gegeven zij f(0) = f(1)= 0 en f © - diff. 
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25e 


| 


Ads 


(n) 


Laat zien, dat de reeks x convergeert in de supnorm naar een 


funktie Xx . 
oo 


Geef een schatting van max xn) _ ln dus van de convergentie- 
nn tel 0,11 
snelheid. | 
Laat m.b.v. de identiteit 
oo : il | 
2 & sin(krt). { Sf sin(kmt) f(T) dr } = Et). 
k=i En: 


zien, dat x inderdaad de gevraagde oplossing van de vgl. + 
randvoorw. 1s- 


Zij s = kâ. 7? met k, € IN. Aán weike voorwaarde moet £ voldoen 


opdat b en c "uitgevoerd" kunnen worden ? 


Beschouw het probleem 


JE + A f(t) x= 0 


[x(0) = Xx(1) = 0 
met f(t) > 0 en f& C° [o,1l. 
Dit probleem heeft uiteraard altijd als triviale oplossing x = 0. 
Als de vgl. + randvoorw. een niet triviale oplossing Xx = vtt) 
bezit voor de waarde à = Aj’ dan heet À … een eigenwaarde en Vn 
de bijbehorende eigenfunktie. | 
Gegeven zij, dat een vgl. + randvoorw. zoals boven een aftelbare 
reeks reêle eigenwaarden Ek SA, S Agee. beate net Än > 


voor n > @, 


‚ Bereken de eigenwaarden en bijbehorende eigenfunkties als f = 1. 


We zullen proberen een benaderinsprocédé van de eigenwaarden in 
het algemene geval op te zetten. 
Geef daartoe een variationele formulering van het probleem en 


tevens het stelsel vergelijkingen waaraan de optimale 








Ln 


Een | | | 
Het da sk = 13..-s n } moet voldoen bij het Kitd-Galerkin 
A 


procédé met als basis een stelsel { 6, 5 X 7 be daasn Heen al. 


(nò) e 
n-de iterand x = Heg dj 


ALI 


bs waarbij &, ven de randvoorwaarden 
voldoet … 


Dit stelsel heeft steeds de triviale oplossing a, = «…… = an * 0: 


Voor speciale waarden van \ echter heeft dit stelsel ook niet 


triviale oplossingen . 


(n) 


Laat zien dat deze waarden Aj volgen uit de vgl. 


m3 


1 ok. 


det [ M 


met Mn) en volk n xn = matrices van de volgende gedaante 


MD (i,3) = le 
1 1,j) = S ò. &. dt voor i f j 
ak 
1 
f (ò.)? dt voor Ì = j 
f ï 
Cad se « _ Î 
M, (1,j) = Sf fCT) d; b. dt voor Ì # j 
0 
| 
ij Et) (6,)° dt voor i = j 


Deze vgl. voor À heeft dus hooguit n opd. 
(n) 


Het valt te verwachten, dat Ar voor n 'groot! onder bepaalde 


omstandigheden een benadering vormt van Ag 
Hoe vindt u een vermoedelijke benadering ge van de bijbehorende 


eigenfunktie ? 


a 
em 


Voer dit procédé uit voor n = 1,2 in het geval dat f = 1 en 


k d 


HH 


Vergelijk uw antwoorden met de exacte A,» Wjs en As Wz uit 


onderdeel a. 


iu. 


On Asymptotiek. 
26. Gegeven zij een continue funktie w op (0,1), die voldoet aan 
vw > 0. We beschouwen het gedrag van py voor t Er 


Verifieer de volgende vergelijkingen 





ae ak = ON S= HD 
ow) 


| b. 1 = OW) > Wb 


01) > exp (-WÌ) = OCH) 


® 
Ee 
Ä 


27. We beschouwen weer het gedrag voor t * 0 van continue funkties 
op (0,1). Definieer het ordesymbool Os als volgt 

| Pz Ogl0) * {Wp = OCH) Één ò = OCH)} 

| a. Laat zien dat 


D= 0g(ô) * JA,B,D > 0 zodat op (0,6): A.lel < vl SB. lol 


Toon verder aan, dat 
1 


14 Os CW 


1 ER wi + Os 


b. p= OC) > (1 WIT 


Sm 2 


) 


OW) at Pv) 


28. De beschouwde situatie situatie zij dezelfde als in opg. 2] 





a. Bewijs, dat de ordening gedefinieerd door het symbool O niet 
totaal is, d.w.z. vind 2 funkties p en & waarvoor beide 
uitspraken : Wp = Old) en & = Ol) ongeldig zijn. 


b. Bewijs, dat het symbool Og een equivalentie relatie definieert. 





29. Gegeven zij een complexe funktie f : C * Cl welke analytisch is 


in een omgeving van de oorsprong. De machtreeksontwikkeling 


ke o „(n) 
f(z) = 2 E_l:,n heeft dus een convergentiestraal p > 0. 


n= 


We zijn geïnteresseerd in het gedrag voor z * 0. 
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30. 


Se 


32. 


Ab: 
Bewijs 
a. de rij {25} 1S asymptotisch geordend 
| En EE B PC zn | 
b. de “eindige Taylorontwikkeling" EF == X end belt 
n=0 ’ | 
een asymptotische ontwikkeling van f voor. 


ö Te Es ocz), 
m 


Ais £M* (0) # g,aan £ - de 


(m+1) 


als f (0) = 0, dan zelfs f - en = o(zAt tj, 


de EES DE HEP, 


Bereken een asymptotische ontwikkeling voor t > 
naar machten van 1/t voor de funktie 
t 
Ò = [U + t.z°) Ee dz met s > 1. 
Sal 


Ga na, dat in dit geval de asymptotische ontwikkeling puntsgewijs 


convergent is en zelfs convergeer naar $. 


Definieer xy en d, als de negatieve, resp. positieve wortel 


van de vierkantsvergelijkirg 


en me Le ta den te 
Laat zien, dat voor t > : X4 = ti + 0(1/+°) 
X) = t + O(1/t) 


Beschouw de transcendente vergelijking 


e“ = tx se Hre Bs 


a. Laat zien aan de hand van een tekening, dat de vergelijking 


twee reele wortels x. en Xj heeft welke voor t > @ voldoen 


1 


aan resp. X4 7 0% Xj > oo , 


b. Bewijs, dat : (t-1) - /TE-D7-7 S x, S H(t-1) - SCEE) 


voor t groot genoeg. 














33. 


34. Beschouw Ct) 





16. 


Das : &, = (t-1)7Î + ot” *) voor t > 

c. Stel x, = log t + y (waarom?). Leid een vergelijking voor 
y af. | 

8 em log lok t : , 

Toon aan, dat : xj log t:. {1 + O0 ( lee )} voor t + 
Let wel : de fout in de schatting X =hen(t) gaat voor 5 
dus absoluut genomen > »@ „ doch relatief genomen > 0. 

d. Toon aan ‚ dat : 5 = log t + log log t + (EEE). 
Merk op dat in deze benadering de absolute je 0 gaat 


voor t > ! 


In het diktaat werd de volgende asymptotische ontwikkeling 


voor t > © afgeleid 


oo m e ee 
fe de, == r-DTK-ijt.t kk 4 Rr 
0 5 k=1 


-(m+1) (*) 


waarbij R_ voldoet aan : IRIS m!.t 
Bereken bij vaste t bij benadering het optimaal aantal termen 
u(t), dat men in de ontwikkeling mee moet nemen opdat Rt) 
geminimaliseerd wordt en geef aan hoe groot en wordt bij dit 
optimaal aantal termen in de ontwikkeling. 


Hint : 1. gebruik de bovengrens voor R, zoals gegeven in (*) 


en benader hierin ms met de formule van Stirling : 
mt Am. ml. 
2. Beschouw |Fmtal/, 
Air 
oo 
=f et Ade dx; t > 0 waarbij f een oo-diff. funktie — 
is op [ 0,9) met Erens | nn 
a. Laat zien, dat de integraal, die Ct) definieert, bestaat. 
b. Vind een asymptotische benadering de) voor ò(t) als t > eo zodanig 


dat de relatieve fout out” ?) is. 
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36. 


der | 


Een stochastische variabele x die waarden aanneemt in R heet star- 
daard normaal verdeeld als de kansdichtheid gegeven wordt door: 
2 
4 =_X 15 
p(x) Jr Oe 


De kans dat x > t is dus: 
| | oo 2 | 
p{x > t} = Te S e”*“/2ax 
E | 


Bereken een asymptotische benadering voor t * ee van deze kans met 


een relatieve fout oct”). 


Beschouw een 9e-diff. funktie f van [0,1] in R. 


a. Laat zien dat f een Fourier ontwikkeling bezit van de vorm: 


oo 
_ EE fa … „èt in X 
aad sE A + 2 XE Re (a, e } 


/ | n=1 
1 | 
met a == f f(x) dx 
o | 
0 | 
1 2m in xXx 
a z=f f(x) … e | e dx 


b. Geef een asymptotische ontwikkeling van a({n) voor n + oo 


ec. Onder welke voorwaarde(n) geldt: 


an) On”) voor n > 
met k CN ? 


d. Bewijs de volgende schatting : 
| (2) an f en 
£ SC [0,1] en f(0) = £(1) » 


bh ee ON bte 
max [f(x) - {a + 2. E Rela, . e@T in Aj < 
0,1] A n 
Ee nk ed A 4 
met Cz — . {I£'COL + |E! + max |£"GO]}. 
| 2 | EE: [ IN 


| 
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le deeltentamen Toegepaste Analyse datum 2 november ; tijd 10-13 uur 


Opgave 1. 


Beschouw de volgende differentiaalvergelijking met beginvoor- 


waarde 

pn dx _ ‚-1 

(1) ae Er W 
(11) A, = x,E R 
met w E R. 


a. Voor welke punten Í x,‚t, | is er volgens een u bekende existentie 
en eenduidigheidsstelling een interval t, <t, < t,, waarop een 


unieke oplossing bestaat van (1) + (ii). 


Bereken de oplossing Xx Ct) van vergelijking (1) welke voldoet 


D. 
aan x, (1) = a ER. ĳ 
Onderzoek het gedrag van deze oplossing voor t + 0. 
Bereken met name lim Xx, Ct) als deze limiet bestaat. 
t+0 | 
c. Ga voor elke B E R na òf er een oplossing van vergelijking 
(1) bestaat welke tevens voldoet aan x(0)= 8 en zo ja, 
hoeveel. 
Opgave 2 


We zoeken een funktie u, welke de volgende funktionaal J mini- 


maliseert : 


_ En du,z 
J(u) = ee, GE) . dt 


en ek 


en tevens voldoet aan de randvoorwaarden 


u(0) = 0 5 ul1) = 1 





| 
t 








ee ne en WE ende RR EE EE EE A LT 


Stel daartoe u = u 


[7 


Laat zien door beschouwing van de bijbehorende Eulervergelijking 
dat er in co, hooguit één u, bestaat welke J minimaliseert 
en aan de randvoorwaarden voldoet. 


Bereken deze kandidaat us 


Laat zien, dat er geen uu. in cr o, 11 is waarvoor. 


1 
Ju) < Jlu,). 


+ n en schat 


1 
1 | 
| _ 5 zl duz dn 


Opgave 3 


Gegeven Z1j 
hi zá ian? 
I(t) = on f (1+2%) 1 e ztz « dz 
… OO 


We zijn geïnteresseerd in het gedrag voor t > @ van I(t). 


/ 


Bereken een asymptotische ontwikkeling van I(t) voor t * @ 


_ 
naar machten van t . 


Î in een Taylorreeks rond z z 0 en 


Ontwikkel daartoe (1+22) 
gebruik de gegevens onderaan deze opgave. 

Noem de benadering die u krijgt door de termen tot en 
met 0(zg) van deze Taylorreeks mee te nemen : 1 (t). 


Laat zien dat er sprake is van een divergente ontwikkeling voor 


vaste t. 


Toon aan, dat 


2n+2 1 


oo 
| _ Alt _AAn+1 Z Jtz? 


Vind een schatting van de gedaante 


(nti) 


(eh HEE TCS Es R_ Ct) 








waarin u zelf Cc heeft berekend. 


Ga bij gegeven t, > 0 na voor welke waarde van n 


minimaal wordt. 


ER Cede 


Beschouw daartoe KR, 


U mag gebruiken, dat 


90 tel hd - En 2 
S „2n, e 24 , dz = (2n-1) Í ze 2 e dz". d 


(2n-1).(2n-3). .….…. … 1e 


KaG) 


Z 


/2n voor n > 1. 


